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Ideas of Statistical Physics are very relevant for cosmic structures especially considering 
that the field is undergoing a period of exceptional development with many new data 
appearing on a monthly basis. In the past years we have focused mostly on galaxy 
distributions and their statistical properties. This led to an interesting debate which will 
be resolved by the next generation of data in a couple of years. In addition here we discuss 
the statistical properties of the fiuctuations of the cosmic microwave background which are 
small in amplitude but complex in structure. We finally discuss the connection between 
these observations and the Harrison-Zeldovich spectrum and its further implications on 
the theories of structure formation and the and the cosmological A^-body simulations. 



1. Introduction 

From the point of view of Statistical Physics cosmic structures give rise to a puzzling 
and apparently contraddictory situation: on one hand the Cosmic Microwave Background 
Radiation (hereafter CMBR) is extremely isotropic with temperature fiuctuations of the 
order 5T/T ~ 10~^ 11)131 • other hand matter distribution is extremely clumpy: 

there are galaxies, clusters, filaments and large voids up to huge scales [^p!^ . Fluctuations 
with a large density contrast 6p/ p ^ 1 appear to characterize the matter field over a broad 
range of scales. The link between the two observations is, of course, not direct. The first 
one refers to radiation which is essentially unchanged from the time of the electro-magnetic 
decoupling. This is when the universe was about 1000 times younger than now. Galaxies 
instead correspond to the present (visible) matter field. Another important difference is 
that the CMBR corresponds to an angular distribution, while the galaxy distribution is 
three-dimensional, in view of the red-shift {z) measurement. Before the decoupling time. 



matter and radiation are coupled by the Sachs- Wo If effect |12| which describes the fact 
that if radiation escapes from a gravitational potential well it is red-shifted by an amount 
linear with the potential. 

This implies that, at z = 1000 (Universe 1000 times younger than now) 6p/ p ~ 6T/T ~ 
10~^. The question is then if gravitational instability can lead to the present large struc- 
tures, starting from Sp/p ~ 10~^ at z = 1000. The answer is no, in fact one could get 
basically an amplification of a factor 10^ by gravity, in the linear regime, leading to a 
present value of the matter fiuctuations Sp/p ~ 10~^, in sharp contrast with the large 
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structures observed. This is one of the main reason to introduce a vast amount of non- 
baryonic dark matter matter (it interacts very weakly with radiation) which is supposed 
to have, at the decouphng time, fluctuations much larger than 10^^. In view of this and 
other puzzling observations the current interpretation is that the matter field consists in 
about 99% of the total matter and the baryonic matter is just about 1%. 

From this discussion it is clear that the understanding of the nature of the various 
cosmic structures and their dynamical development is a central point of the field which 
has implications in essentially all other properties. Many new data have been generated 
in the past few years and many more are expected in the near future. So the field is 
quickly evolving from a situation with many speculations and few data to a new one, in 
which the new abundance of data permits a test of long standing conjectures. This is 
therefore an extremely interesting moment in which ideas from Statistical Physics can 
play a crucial role. This is true grasping the properties of structures from the data but 
also in the theories of structures formation and in the A^-body simulations. 

In the past we have mostly considered the properties of the galaxy distributions 
Here we start instead from an analysis of the fluctuations corresponding to the microwave 
background radiation, and then we try to connect them to the galaxy properties. We out- 
line the main open problems from a new perspective and we consider some consequences 
also on the theories and cosmological A^-body simulations @J^. 

2. Small fluctuations and Poisson distribution 

The standard approach in the field is to assume that the matter field p(r) is always 
characterized by a well defined average density with small amplitude fluctuations. In our 
opinion this may be acceptable in relation to the microwave radiation, but not for the 
galaxy distribution. We will come back on this point later. 

The basic assumption is to have 

p(f) = po + Spir) with 5p{r) < po , (1) 

where po is the average density. This is then analyzed in terms of the correlation function 
of the fluctuations (positive and negative) 

^ _ (p(ro)p(fo + r)) _ ^ {5p{ro)5p{fo + r}) 
Po Po 
the normalized mass variance in spheres of radius R 

^ ^ {M{R)f ' ^ ^ 

(where M{R) is the mass in a given sphere of radius R) and the power spectrum of the 
normalized density fluctuations which are supposed to have small amplitude with respect 
to the average density 

P{k) = {Mk)\') = ^ / rfVe-n(r) , (4) 



where 6p{k) is the Fourier transform of the normalized density contrast field 6p{r)/po. 
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It is very instructive to consider various properties and in particular the fluctuations of 
the gravitational potential corresponding to a Poisson distribution of mass points. The 
Poisson distribution is completely random and uncorrelated. In a sphere of radius R we 
have an average number of points given by 

{N{R)) = poV PoR\ (5) 

where V ^ R^ is the volume of the sphere in three dimensions. The typical number 
fluctuation in such a sphere is given by 

5N{R) = {N{R))^ ^ R^ . (6) 

With respect to our functions we have that in the thermodynamic limit of an infinite 
distribution of points 

^(r) = 6{f) ; P{k) r^k'^ ^ const. > . (7) 

However in a finite sample, for the finite size fluctuations due to Eq.(|^), for a generic 
statistical estimator of ^(r) we have 



6p 




6N 


Po 




{N) 



{Ny^ ~ r-t (8) 



for the absolute value of ^(r), where 6p is the typical fluctuation of the conditional density 
from an occupied point Q For P{k) the finite value of N induces some noise on top of the 
constant value of Eq.(0). 

Now comes the key-point because we can estimate the fluctuations of the gravitational 
potential S(j){R) function of the scale R: 

5<P{R) ~ ~ R^ . (9) 

R 

So we have diverging fluctuations of the potential originated by the finite size fluctua- 
tions of the Poisson distribution Eq.(^. This is a very important point which has deep 
implications with respect to the Harrison- Zeldovich spectrum (see below). However, if 
we estimate in the same way the fluctuations of the gravitational force F{R) we have the 
surprising result that they are instead bounded 

\F{R)\^\V5(t){R)\^ R'^ . (10) 

This result, that the gravitational force is well defined in a Poisson distribution of points, 
is well known since the seminal work of Chandrasekhar [^. The extension of this result 
to the fractal case is rather complex, and the interested reader can find a discussion in 

i- 

These properties of the gravitational field of a random distribution of of massive points 
are usually not very appreciated but they already pose, in our opinion, extremely inter- 
esting and challenging problems as we are going to see in the following. 



^Note that there are some special statistical estimators of ^(r), for the particular case of a Poissonian 
distribution only, giving |^(r)| ~ r~^. 
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3. The Harrison-Zeldovich (HZ) spectrum 

We consider now the standard derivation of the HZ spectrum of primeval fluctuations 
which has a central role in the field. The idea is to consider the implications of the 
CMBR observations as given for example by the COBE experiment . In this experiment 
the temperature fluctuations of the black body spectrum are measured at large angular 
separation ranging from 7 to 90 degrees. The observation is that the value 5T /T ~ 10~^ 
is about the same for small and large angles. This means that this fluctuation does not 



grow with the angle. The standard conjecture (i.e. [|T2|) is that one can go from angles 
to distances and argue that 

5T{e) 5(I){R) r-. const . (11) 

Namely, if the temperature fluctuation is angle independent, then the gravitational po- 
tential fluctuation which, via the Sachs- Wolf effect led to the temperature fluctuation, 
should be independent on the scale. This is sometimes defined as a conditional "scale 
invariance", but instead is just that the potential fluctuations are simply constant at 
various scales We can immediately see that a Poisson distribution does not satisfy 
this condition because the fluctuations of the gravitational potential grow with scale as 
in Eq.(^. It is easy to see that, in order to satisfy the Eq. (|Tll) , one should have 

5N{R) ~ (iV(i?))^ , (12) 

which implies a more regular distribution {super-homogeneous) with respect to a random 
Poisson one (see Fig. |1|). With respect to the correlation function and the power spectrum, 
it can be easily shown that at large scales ^ 

iir) ~ -r"^ (13) 
P{k) ~ k. (14) 

Note the negative sign for the large scale behavior of ^(r): this means anti-correlations 
at large scales (see Fig. |^). 

This discussion leads to a series of questions about the implications of the density fluc- 
tuations on the gravitational potential which, in our opinion, would require some deeper 
clarifications: (i) a Poisson distribution has a perfectly well defined average density at 
large scales po = limy^oo ^(^)/^- Consequently, one could suppose that at large scale 
one can use the value po iii general relativity and describe it in terms of the Friedman 
metrics. However, given our discussion and Eq.(D, this should not be permitted because 

1 /2 

the finite size fiuctuations 5N ~ (A^) ' induce diverging gravitational potential fiuctua- 
tions even at large scales where 5N <^ (N). If, on the other hand, one would compare 
forces and not potentials, then the perturbation would be negligible. So the question to 
clarify is whether a Poisson distribution of mass points is compatible with a Friedman 
cosmology at large scales, (ii) The same reasoning holds for the derivation of the HZ 
spectrum as necessarily implied by the observed properties of the CMBR. Note that the 
HZ spectrum is so central to the field that this and the previous points should either 
be reinforced on disproven. (iii) Other points concern the relations between angular and 
space correlations. We have learned that for galaxies the two properties can be very dif- 
ferent due to the projection effects of the angular case. Finally the Sachs- Wolf effect for 
non trivial density fluctuations should be considered. 
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Figure 1. Left panel: Glass-like distribution of points obtained by the method of the 
shuffled lattice . The idea is to decrease the fluctuations below to the level of the Poisson 
random distribution. This achieved by imposing a maximal value for the nearest neighbors 
distance. In this case the fluctuation is about SN{R) ~ {N{R))^^^ and it is a qualitatively 
good representation of the HZ spectrum which in this respect can be considered super- 
homogeneous because it is more homogeneous than a Poisson distribution. Right panel: 
Example of a Poisson distribution. In this case the behavior of the typical fluctuation of 
the number of points in a ball of radius R goes like 6N{R) ~ {N{R))^^^. The distribution 
is homogeneous at large scale with a well deflned average density. 



4. Galaxy correlations updated 

For galaxy correlations all the above concepts of small fluctuations theory etc. are 
certainly not appropriated. One observes fractal clustering with large scale structures and 
voids for which 6p/pQ^ 1 or even worse po is not well deflned in the available samples. 
This is the fleld in which we have worked extensively and proposed a radically new and 
more general way to interprete the puzzling observations and to deflne the theoretical 
perspective. One of the consequence of our work is that the so called characteristic length 
of galaxies tq ~ 5Mpc/h |jl2| is an artefact of the a-priori assumption of homogeneity 
and, in our opinion, correlations extend up to the present limits of observations P, p!^ , p 



In any case there is, flnally, agreement on the fact that at least to some small scale, 
galaxy correlations are fractal. A common misunderstanding is however to think that if 
there would be a crossover, then the small fluctuations theory could describe the fractal 
properties because they extend only on a flnite range of scale (e.g. [15 1). This is clearly 



incorrect because, independently of the presence of the crossover, a theory based on small 
fluctuations with Sp/po -C 1 can never be appropriated to describe fractal clustering. 
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Figure 2. Left panel: Power spectrum corresponding to the HZ spectrum (solid line). 
Note that the behavior at large scale (small k) P{k) ~ k (dotted line) implies less fluc- 
tuations than the Poisson behavior P{k) ~ ~ const. > (dashed line). Right panel: 
Modulus of the correlation function ^ (r) (solid line) corresponding to the previous power 
spectrum. Note that at large scales the function is negative (the reference dotted line has 
a slope r~'^). Note that the units are arbitrary. 



Apart from this conceptual point, which is anyhow valid, the question of the extension 
of the fractal properties should be clarified by next catalogues. In the next few years the 
new catalogs like 2dF and SDSS will dramaticaly improve our knowledge of the large scale 
distribution of galaxies. The ongoing debate in the field P,|l3|,|l5|,|lT| will be then clarified 
by the analysis of these new data. 

It is clear therefore that, at least for the properties of galaxy clustering, the formalism 
of the small fluctuation theory is not appropriate at any stage from the data analysis to 
the formulation of a theory of structure formation. 



5. Conclusion and perspectives 

We have presented a brief overview of the relevance of Statistical Physics for cosmic 
structures. Many new data are appearing for the CMBR and for the galaxy distribution, 
and modern Statistical Physics is the natural framework to identify the relevant structures 
and defining the corresponding theoretical models. 

In this respect the situation appears quite different for the very small fluctuations in 
the CMBR and the very large fractal fluctuations observed for galaxies. In the first case 
the traditional theory of small fluctuations could be appropriated, but it is certainly not 
appropriated for galaxies. Whatever clustering is observed this has scaling and fractal 
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properties and, even if there would be a crossover to an homogeneous distribution (a fact 
which is still far from shown) the fractal clustering could never be properly described by 
a small amplitude fluctuations theory. 

A different question concerns the cosmological model able to describe a situation where 
matter has fractal properties. In Ref. |T0[| the interested reader can find an attempt to 
construct a Friedman model with zero density, where the fractal can be treated as per- 
turbation to the background homogeneous metric. 

It is therefore unavoidable to generalize the approach to include the possibility of com- 
plex structures and large fluctuations at alia levels, from the data analysis to the the- 
oretical methods and the A^-body simulations This should lead to a new general 
perspective of the field, strongly based on a cross- disciplinary approach, which certainly 
will lead to novel exciting developments. 
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